Dynamical coupling and separation of multiple degrees of freedom in a photoexcited 

double-exchange system 
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We present a theory of ultrafast photo-induced dynamics in a spin-charge coupled system, motivated by pump- 
probe experiments in perovskite manganites. A microscopic picture for multiple dynamics in spin and charge 
degrees is focused on. Real-time simulations are carried out by two complimentary methods. Our calculation 
demonstrates that electron motion governs a short-time scale where charge and spin dynamics are combined 
strongly, while, in a long-time scale controlled by spin relaxation, charge sector does not follow remarkable 
change in spin sector. Present results are in contrast to a conventional double-exchange picture in equilibrium 
states. 
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Ultrafast photo-control of electronic and magnetic struc- 
tures has attracted much attention for a long time from view 
points of fundamental physics and technological application. 
In correlated electron systems with multiple degrees of free- 
dom, i.e. charge, spin, orbital and lattice, a stable equilibrium 
phase is determined by a subtle balance of several interac- 
tions between them IH. In a barely stable state at vicinity of 
phase boundary, a gigantic change in electronic structure is 
triggered by an optical pump pulse. Recently developed sev- 
eral time-resolved experiments enable us to access directly to 
the photo-dynamics of the multiple degrees of freedom 

Exotic equilibrium phenomena in correlated electron sys- 
tems, such as, high-Tc superconductivity, colossal magnetore- 
sistance (CMR), multiferroics, and others, are often attributed 
to strong correlation between multiple degrees of freedom. It 
is not quite trivial whether this naive strong-coupling picture 
between multiple degrees is applicable to the photo-induced 
dynamics or not. One example is seen in perovskite mangan- 
ites ll, 3, S 0] . A key issue is a coupling between spin and 
charge. The CMR effects are resulted from a subtle balance of 
the charge ordered (CO) insulating phase associated with the 
antiferromagnetic (AFM) order and the ferromagnetic metal- 
lic (FM) phase, and are addressed in strongcoupling between 
localized spins and conductive electrons yj]. By irradiation 
of a femto-second laser pulse in the CO insulating phase near 
the boundary, a transient metallic state expected from a re- 
flectivity change appears within lOOfs. While a macroscopic 
magnetization suggested from the magneto-o ptic al Kerr rota- 
tion increases with the time of few ps IITilsllgillOl 11]. 

These results imply that changes in spin and charge sec- 
tors, which are recognized to be strongly correlated in man- 
ganites, are not always accompanied with each other in photo- 
excited state. This dynamical coupling and separation of mul- 
tiple degrees of freedom is not only a clue for an overall un- 
derstanding of the optical-control electronic phases in corre- 
lated system. This also provides a key to uncover exotic phe- 
nomena in equilibrium state; roles of the multiple degrees are 
able to be identified separately. In spite of the progressive re- 
searches, however, a microscopic picture of multiple dynam- 



ics in photo-excited states has not been clarified yet. 

In this Letter, motivated from the pump-probe experiments 
in manganites, we present a theory of the photo-induced dy- 
namics in a spin-charge coupled system described by the 
double-exchange (DE) model. We focus on dynamical cou- 
pling and separation between spin and charge. To attack this 
issue, two complimentary methods, the exact-diagonalization 
(ED) and the Hartree-Fock (HE) methods, are adopted. Cal- 
culations demonstrate that electron motion governs the short- 
time scale where charge and spin sectors change coopera- 
tively. On the other hand, in the long-time scale controlled 
by spin relaxation, charge sector does not follow change in 
the long-range spin correlation. 

Let us introduce the extended DE model. 



DE = -at [cl^Cja + H.c.j -JnY^Si-Si 
{'])" ' 
+ UY,n^nii+VY,ninj+JsY,Si-Sj, (1) 

m (ij) 

where c,a is the annihilation operator for the conduction elec- 
tron at site ; with spin fl(=t, j), and S,- is the operator for 
the localized spin. We introduce the number operator «, (= 
laW/a = L,4,Ci«) and the spin one s,- = 5 Lfliq„c7„foC,fo with 
the Pauli matrices a for the conduction electrons. The first 
and second terms describe the electron transfer, at, and the 
Hund coupling, Jh, respectively. The on-site Coulomb re- 
pulsion U, the nearest-neighbor (NN) one V , and the AEM 
superexchange interaction Js between the localized spins are 
also taken into account. All energy and time parameters are 
given as a unit of t and f^^ respectively. A magnitude of 
the transfer integral is changed by changing a from one. Eor 
a typical value of t in manganite, a unit of time, f^' ,corre- 
sponds to O.Sfs'-^lfs. 

To analyze this model, two complimentary methods, the 
ED and HE ones, are utilized. In the ED method, the time- 
dependence of the wave function and the excitation spectra 
are able to be obtained exactly. On the other hand, in the HE 
one, simulations in a large system size and long time distance 
are possible, and the spin relaxations are taken into account. 
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We utilize, in the ED method, the time-dependent and dy- 
namical density-matrix renormalization group (DMRG) and 
Lanczos methods 11211 . One-dimensional clusters of system 
size Ni ~ L{< 13) with the open-boundary condition are 
adopted. Electron number is Neu = (^+ l)/2 corresponding 
to the 1/4 filling. For simplicity, amplitude of the localized 
spin is set to be 1 /2. We take a damped oscillator form for the 
vector potential of the pump photon as Ap,„„p = Aoe^"**^^^l^l 
at time T with frequency Oo and a damping constant 70. A cen- 
ter of the wave packet is defined by T = 0. The wave function 
for the one-photon absorbed state at time t(> Jq^) is derived 



by the first order perturbation with respect to A pump 11 1311 : 



7o 



(coo-^DE + Eor + Y^ 



m, (2) 



where j = iatY,{ii)aWia^ia—H.c.) is the cuiTent operator, |0) 
is the wave function before pumping, £0 = (0|.^^£|0), and 
,_A^ is a normalization factor. Transient excitation spectra are 
calculated by the linear-response theory. For example, the 
one-particle excitation spectra are given by a sum of the elec- 
tron and hole parts, A(q', (o) = A^^^{q, (o) +A'"''''{q, co), where 
the first term is obtained as iu 



A"'' (a, CO) 



1 



(o-JfoE+E + iy 



Cqa)- (3) 



We introduce the Fourier transform (FT) of the operator, c^a, 
the energy E = (J^de), and a damping factor 7. A bracket 
(• • • ) implies the expectation value with respect to |0) in the 
case before pumping, and that with |1'(t)) after pumping. The 
pseudo-momentum in the open boundary condition is defined 
as Og = y^2/{L + 2)Y,i sin{qri)Oi where (9,- ((9^^) is a physical 
quantity in the real (momentum) space. 

In the HE method, the mean-field (ME) decoupling 
is introduced in the many -body terms in J^foE, such as 

const, where (n,|) is a site-dependent ME. The two- 
dimensional Nl — L X L (L = 6 — 10) site cluster with the 
periodic-boundary condition is adopted. The electron num- 
ber is Neie = L^/2. The localized spins are treated as classical 
vectors with amplitude of 3/2. The pump-photon irradiation 
is simulated as an electronic excitation from the highest occu- 
pied HE level to the lowest unoccupied one without changing 
spin. The time evolution of the electronic wave function is 
obtained by the time-dependent HE scheme 1 15 , 3 ; 



^<'^(t + At) = Texp 



t+At 



0«(t), (4) 



where 0(t)('' is the /-the HE wave function, and -^^1^/ {t') is 
the HE Hamiltonian for Eq. ([T]). Dynamics of the localized 
spins is described by the Bloch-type equation: 



drSi — Hi X Si + d^Silrelax, 



(5) 



where H, is the ME at site / defined by H, = ///(s,) + 
JsHs^i+S- The last term in Eq. Q represents the spin re- 
laxation and dephasing which will be explained later. 
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FIG. 1: (color online) (a) One-particle excitation spectra before 
pumping, (b) those at rt = 3, and (c) those at Tt = 10. Momenta 
are changed from ;r to for the data from the top to the bottom, (d) 
Time-dependence of optical absorption spectra. 



Eirst, we show the results obtained by the ED method. The 
parameter values are chosen to be f/ = 10, /// = 8, V = 5, 
Js = 0.4, too = 3.72 and 70 = 0.4. We take a ~ \ except 
for the results in the inset of Eig. [3] Somewhat larger en- 
ergy parameter values than the realistic ones for manganites 
are used to reproduce the AEM/CO insulating ground-state in 
one-dimension. Thus, our analyses in the ED method are re- 
stricted to qualitative properties in the photo-induced dynam- 
ics. Eirst we focus on the charge dynamics. An insulating 
gap and a tendency of the Brillouin-zone doubling shown in 
A((i, (a) [see Eig.[TJa)] imply an alternating charge alignment. 
Just after photon pumping [see Eig.[TJb)], a photo-carrier band 
appears inside of the CO gap in A(^, (o). Large amounts of the 
spectral intensity in the upper and lower bands still remain. 
These results imply that the long-range charge alignment is 
collapsed, but a short-range correlation survives. It is worth 
to note that the band width of the in-gap band is broaden with 
time evolution [see Eigs.fllb) and[Tlc)]. The optical absorp- 
tion spectra a{o}) = -{nNLy^lm{j{co - ,J%e + E + iy)^^ j) 
are presented in Eig. [lid). After photon pumping, spectral 
weight appears inside of the optical gap, and grows up with 
increasing in time. We confirm, through the analyses in a 
small cluster system, that the lowest component of the in-gap 
spectral weight corresponds to the Drude component in the 
thermodynamic limit. These time evolutions in A{q,(o) and 
a{(o) are not observed in both the spin-less V-t and Hubbard 
models. That is, the localized spins and its coupling with con- 
duction electrons play a central role on these phenomena. 

Change in the spin sector is monitored by the static spin cor- 
relation function, ^(q) = (S_q • Sq), where S^, is ET of S, [see 
Fig. 121. Before pumping, the spin structure is ferrimagnetic 
one where the localized spins are aligned antiferromagneti- 
cally and the total spin is S'"' = 5 + This is seen in large 
intensity of S{q ^ n). After pumping, the large AE correla- 
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FIG. 2: (color online) Spin-correlation functions at several times. 

tion is rapidly suppressed, the momentum dependence in S{q) 
is almost smeared out, and the system is nearly paramagnetic. 

Time evolutions of both the conduction electron and local- 
ized spins are summarized in Fig. [3] We plot the second mo- 
ment W of A {q,(o) reflecting the band width of the in-gap 
state, the integrated spectral weight D inside of the optical 
gap in a{co), and the NN spin correlation Ks for the local- 
ized spins. These are defined by W = f^" ^j^A(g', co)(a) — 
(Oc)^d(0 with the center of the in-gap band (Oc and the up- 
per (lower) band edge couicoL), D = f'^f^ a((o)dco with the 
upper (lower) edge of the in-gap component co^{co[), and 
Ks =N,j^ (S/ • Sj) with the number of NN bonds Nb- An 
almost identical time dependence is shown in the three curves; 
they increase linearly after pumping and are saturated around 
T = IQ/t. The characteristic time scale where Kg shows a 
shoulder is denoted as T^. In the inset of Fig. [3] Kg's for differ- 
ent values of the electron transfer are presented. With increas- 
ing a, a slope of the curve in the region of T < Tj increases, 
and the AFM correlation collapses rapidly. When Kg's are 
plotted as functions of rat, the slopes for several values of 
a are almost identical. We conclude that photo dynamics 
in the charge and short-range spin sectors are strongly corre- 
lated with each other, and are mainly governed by the electron 
transfer These results are interpreted as follows; motion of 
the conduction electrons, excited directly by photon, destroys 
the charge and AFM orders. The collapse of the AFM spin 
alignment further promotes motion of the conduction electron, 
and finally the two systems steady down cooperatively. 

In the calculations by the ED method explained so far, the 
spin-angular momentum is conserved. However, it is ex- 
pected that, in the simulation for longer time scale, the spin 
relaxation leads the system to the lower excited state on the 
energy surface. Microscopic mechanisms for breaking of 
the spin conservation in manganites are not clear yet. In 
photo-excited states where the Mn 3de^ orbitals are mainly 
concerned, the orbital angular momentum, L, is basically 
quenched. This is in contrast to the magnetic semiconductors, 
such as (Ga,Mn)As where a large spin-orbit coupling provides 
a spin relaxation UtIi . Possible relaxation mechanisms are 



FIG. 3: (color online) Width of the in-gap band in one-particle excita- 
tion spectra, W, spectral weight inside of the optical gap, D, and NN 
spin correlation function. Kg. Data are subtracted by values at T = 0, 
and are normalized by differences between the minimum and maxi- 
mum values. We take (tOj/,©^) = (6.5,-1) and (&)[,, a)[) = (4,0.5). 
The inset shows A^^ for several transfer integrals. 

mixing of the on-site eg and t2g orbitals by the off-diagonal 
components of L, and that induced by lattice distortions with 
the Tiu symmetry. A mixing of the inter-site and t2g or- 
bitals by the GdFeOs-type distortion is another candidate. We 
predict that the photo-induced magnetization is promoted by 
excitation of the phonons concerning these lattice distortions. 

We introduce the spin relaxation phenomenologically in the 
HF formulation, and examine roles of the spin conservation 
on time evolution. In Eq. (|5]l, we introduce the parallel and 
perpendicular spin components to the total magnetization S = 
eS, with a unit vector e, as II17I1 S/ = S,-^ + Sqe. The relaxation 
terms are given by d^Si\\ \reiax = -^L{Si\\ -S) and d^Si^\reiax = 
— FrS;^, where and Fj- are the relaxation constants with a 
relation Fz. = 2F7-. A value of f^tir) should be much small 
than that for (Ga,Mn)As where the spin-relaxation is expected 
to occur within lOOfs. We also examine the Landau-Lifshitz- 
Gilbert-type relaxation d-^Silreiax = ^LG^i ^ (^tS,, and confirm 
that results are not sensitive to the types of the relaxation. 

The results by the HF method are presented in Fig.|4] The 
parameter values are [/ = 8, 7// = 6, V = 1 and Js = 0.03. We 
take a = I except for the inset of Fig. EJa). Before pump- 
ing, CO associated with the canted-AFM alignment in the lo- 
calized spins is realized. Two time scales are recognized in 
the NN spin correlations for the localized spins Kg shown in 
Fig. St a); a sharp increase in a short-time scale (~ 200/f) de- 
pending on the transfer integral [see the inset of Fig. Ufa)], 
and a slow increase in a long-time scale depending on F^. 
The short-time scale corresponds to T5 introduced previously. 
The long-time scale termed Tf, is defined as a time when an 
amount of change in the FM spin correlation 5(0) becomes 
1 /e of its total value; Tl = 645 /f and 250/f for F^ = 0.002 
and 0.008, respectively. In the case without relaxation, a Ks 
v.s. T curve is qualitatively similar with that in the ED method 
[see Fig. O. This implies that with regard to the short-time 
dynamics with F^ = 0, a qualitative feature commonly ob- 
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FIG. 4: (color online) (a) NN spin correlation function for several 
values of the relaxation constant F^. The inset shows Kg for the 
several transfer integrals at and r£ = 0.002. (b) Spin correlation 
functions 5(0) for several F^ (bold lines), and S{7t, 7t) forF/, = 0.002 
(dotted line). The dotted arrows represent T/,. (c) Charge correlation 
function N(7t, 7t) and NN charge correlation Ajy. 

served by the two methods is reliable. A quantitative differ- 
ence between the values in the ED and HF calculations 
is mainly attributed to the different photon densities x,,/,; the 
present Xp^ = L^^ = 0.01 is smaller than Xp/, ^ 0. 1 in the ED 
calculation. For small Xp/,, a system takes long time to be set- 
tled down to a stationary state. We checked the above points 
by the HF calculations in one-dimension; a linear increase and 
a saturation are observed in Ks, and Ts increases by increasing 
L from 10 to 30. 

We focus on the results for Fi = 0.002. As shown in 
Fig. Efb), around T5, S{k, k) disappears, but only a subtle 
change is seen in 5(0); the system becomes almost para- 
magnetic. In the long-time scale beyond T5, a remark- 
able change is seen in 5(0) which grows up monotoni- 
cally toward the full-spin polarization, in contrast to Ks and 
S{n,n). Charge dynamics, i.e. the long-range charge cor- 
relation A^(q) = (A«qA«_q) and the NN correlation Kn = 
(2/A^B)^^,y) (A«,An;), are shown in Fig. |4] (c). We define 
A«, = «; — 1/2 and its FT, A«q. Time dependences of the long- 
and short-range charge correlations are similar to that in Kg, 



not in 5(0). These are nearly saturated around and do not 
show a remarkable change beyond T5. Changes in 5(0) by spin 
relaxation are almost irrelevant to the charge sector. This weak 
correlation between the charge correlations and 5(0), corre- 
sponding to magnetization, is not expected from the equilib- 
rium state where these are combined strongly. The key ingre- 
dient in this dynamical separation between the spin and charge 
sectors is spin conservation; magnetization dynamics is gov- 
erned by the spin relaxation rather than the charge motion. 
The present case is in contrast to the photo-induced dynamics 
in electron-phonon systems where conserved quantities in a 
phonon sector does not exist. 

From the calculated results, we provide a whole picture of 
the photo-induced spin and charge dynamics. In the short- 
time scale, charge motions governed by the electron transfer 
destroy the CO and short-range AFM correlations coopera- 
tively. Excess energy given by the pump photon flows from 
the conduction electrons to the localized spins. This state is 
further relaxed through the spin relaxation in long-space and 
long-time scales. Changes in the charge sector is almost satu- 
rated in the short-time scale and do not follow the spin relax- 
ation. This scenario provides a microscopic interpretation for 
the experimentally observed multiple time scales in AR and 
AO. The present dynamical coupling/separation in the mul- 
tiple degrees of freedom is not expected from a conventional 
picture in equilibrium state and are key issues to reveal the 
photo-induced dynamics in correlated system. 
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